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ABSTRACT 


This  thesis  describes  the  development  of  a  new  method 

for  controlling  the  error  in  model  following  control  systems. 

The  treatment  is  for  first  order,  linear  or  nonlinear, 

time  varying  or  time  invariant  systems  with  additive 

(linear)  control.  The  errors  controlled  are  assumed  to 

have  arisen  from  external  disturbances  or  from  differences 

in  the  initial  conditions  of  the  plant  and  the  model. 

The  theory  introduced  is  called  constant  error  model 

following  control.  This  paper  describes  the  theory  as  an 

outgrowth  of  attempts  to  control  a  plant  by  feedback  of  an 

error  between  the  plant  and  model  with  the  error  specified 

« 

to  be  constant.  From  this,  it  is  shown  that  a  model  may  be 
followed  with  arbitrary  error.  The  central  result  is 
that,  given  some  error,  one  can  find  another  model  (control 
model)  which,  if  followed  with  this  arbitrary  error,  will 
guide  the  plant  state  trajectory  back  to  that  of  the  model. 

A 
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I.  INTRODUCTION 


The  theories  presented  in  this  thesis  were  developed 
in  the  course  of  researching  control  techniques  for  use  in 
aircraft  departure  prevention.  The  aim  was  to  apply  model 
following  control  methods  to  the  nonlinear,  time  varying 
and  coupled  dynamic  equations  of  motion  of  an  aircraft  at 
critical  combinations  of  angle  of  attack  and  angular  rates. 
When  it  was  seen  that  a  new  development  in  model  following 
control  had  evolved  from  the  research,  it  became  the  sole 
subject  of  the  thesis. 

Model  following  control  is  concerned  with  techniques 
which  cause  a  physical  plant  to  behave  as  much  like  a  model 
as  possible.  Motyka  [Ref.  1]  accomplishes  this  by  solving 
the  plant  state  equations  for  their  controls  and  then 
substituting  expressions  for  the  model  states  and  state 
rates  into  these  equations.  That  is,  assuming  that  the 
plant  is  defined  by  the  linear,  small  perturbation  constant 
coefficient  differential  equation: 

*P  =  £i  xP  *  *i  UP  (1) 


where 


is  the  state  of  the  plant 
is  the  control 
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and  that  a  model  is  given  by  the  corresponding  differential 
equation: 


x_  *  f,  xm  +  g,  u  (2) 

m  L  m  °2  m 

It  is  desired  that  x  =  x„  and  x  =  x 

m  p  m  p 

The  plant  control  which  achieves  this  can  be  determined  by 
substituting  the  desired  relationships  into  the  plant 
equation: 


x  *  f,  x  +  g,  u  (3) 

m  1  m  &1  p  v  ' 

And  solving  for  u^  .  The  result  is  an  expression  for  the 
plant  inputs  which  make  the  plant  state  equal  to  that  of 
the  model. 

Problems  with  this  method  arise  when  equations  (1) 
and  (2)  are  vector  equations  with  fewer  controls  than  state 
to  be  controlled.  Also,  if  equation  (1),  which  is  a 
mathematical  description  of  a  real,  physical  plant,  fails 
to  describe  that  plant  accurately,  errors  may  be  introduced 
into  the  response  of  the  system.  Finally,  errors  will 
occur  if  the  plant  and  the  model  do  not  have  the  same 
initial  conditions,  or  equivalently,  in  the  presence  of 
external  disturbances. 

The  first  two  of  these  three  problems  are  treated  by 
Moytka  [Ref.  1]  and  Rynaski  [Refs.  2  and  3].  They  are 
beyond  the  scope  of  this  paper.  The  theories  presented 
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herein  apply  specifically  to  the  problem  of  reducing  the 
system  error  in  the  third  case.  The  method,  called  constant 
error  model  following  control,  is  developed  and  demonstrated, 
for  nonlinear,  time  varying,  first  order,  single  input - 
single  output  systems  with  additive  (linear)  control. 
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II.  PROBLEM  DESCRIPTION 


For  the  purposes  of  this  discussion,  it  is  assumed 
that  a  plant  (physical  process)  is  given,  and  that  it  is 
completely  described  by  a  first  order  differential  equation 
with  additive  control  of  the  form 


=  f.  (xn  ,  t)  +  g. 


where  is  the  state  of  the  plant 

t  is  the  variable  time 
Up  is  the  control  input  to  the  plant 
f-^  is  (in  general)  time  varying  and  nonlinear  in  Xp 
g^  is  constant 

It  is  desired  that  the  response  of  the  plant  be  modified  in 
some  way.  This  modified  response  is  completely  described 
by  the  mathematical  model 


'  f 2  (xm  -  tJ  *  8; 


where  x  is  the  state  of  the  model  which  has  the  desired 
m 

response 

us  is  the  control  input  to  the  system  which  incor¬ 
porates  the  model  and  the  plant,  and 


f2  is  (in  general)  time  varying  and  nonlinear  in 

xm  ,  and  f2  f  fj 
g2  is  constant,  g2  /  gx 
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Ill .  SYSTEM  DESCRIPTION 
A  parallel  model  following  system  is  chosen: 


FIGURE  1.  PARALLEL  MODEL  FOLLOWING  CONTROL  SYSTEM 
where  e  is  the  magnitude  of  the  system  error: 

e  *  |  x  -  x  J 
1  m  p 

K  is  the  adaptive  gain  applied  to  e  to  generate 

cL 

UK 

is  the  (additive)  modification  to  the  system 
control,  u 

s 

In  this  system,  the  steady-state  error  (ess)  can  never  be 
zero  unless  the  plant  and  model  have  identical  responses 
to  the  same  input,  since 

ess  =  0  implies  that  u^  =  0  implies  that  Up  *  ug 


ess  =  0  imPlies  that  xmss  "  xpss 


IV.  CONSTANT  ERROR  MODEL  FOLLOWING  CONTROL 


A.  CONSTANT  ERROR 

Without  changing  the  system,  the  error  (e)  is  specified 
to  be  fixed  at  some  arbitrary  value  (e): 


FIGURE  2.  CONSTANT  ERROR  MODEL  FOLLOWING  CONTROL  SYSTEM 

We  wish  to  determine  the  gain,  K*  ,  which  will  insure  that 
this  condition,  once  established,  will  be  maintained  as 
shown  in  Figure  3. 


FIGURE  3.  PLANT  AND  MODEL  STATE  TRAJECTORIES 

The  state  trajectory  of  the  plant  follows  that  of  the  model 
with  constant  difference  e  . 

i 


We  have  defined 


e 


A 


x 


P 


(4) 


so  that  e  =  lx  -  x  *|  =  constant, 

'  m  p  1 

where  x^*  is  the  plant  output  for  e  =  e  .  Similarly, 
denote  the  value  of  K  which  maintains  the  error  constant 
at  e  =  e  as  K*  and  the  resulting  control  as  u*  . 

We  now  have 


e  * 

i  xm  -  X  * 

1  m  p 

i  = 

(xm 

-  Xp»)  sgn  (xm  - 

X  *) 

P 

c  n 

Y  *  = 

Xm 

sgn 

(x  -  X  *)  -  £ 
m  p 

P 

sgn 

(x  -  x  *) 
m  p 

or 

X 

xs 

Ji¬ 

ll 

X 

m 

-  e 

sgn  Cxm  -  xp*) 

(5) 

and 

if 

o 

II 

(x 
v  m 

-  X 

P 

*)  sgn  (xm  -  xp*) 

fat  txm 

-  y>i - 

0 

since  (x  =  x  *)  is 
v  m  p  ' 

assumed 

constant  ) 

so 

x  * 

P 

=  X 

m 

IS 


\ 


The  esult  is  that,  if  the  plant  state  is  exactly  +_  r 

from  the  model  state,  we  may  express  x  *  and  x  *  as 

P  P 

functions  of  x  ,  x  and  +  e  .  The  values  of  x  ,  x 
mm-  m  m 

and  e  are  known  for  any  time  and  system  control  because 
they  have  been  specified. 

Now  note  from  tne  system  description  (Figure  2)  that 

u  *  =  u  +  u  * 
p  s  k 


or 


u,,  *  =  u  *  *  u 
kps 


(7) 


and 


u. 


*  =  v  * 


k  *  E 


(8) 


Equations  (1)  and  (2)  may  be  solved  for  u^*  and  ug 


*  '  (fx  (x  *  .  t) 

u  *  =  _E - ± - 2 - 

P  gx 


(9) 


u 


*  -  f.  (x  ,  t) 

m  2  m 


(10) 


By  substituting  equations  (5)  and  (6)  into  equation  (9) : 


-  fi  [(xm  -  e  sgn  (xm  '  x  *))  ,  t] 


u 


m 


m 


m 


Si 


(ID 


If  U  and  u  are  not  additive  controls  as  assume 
P  S 

and  if  the  state  equations  may  be  solved  for  u  and  u. 
the  results  which  follow  are  still  valid.  P  ' 
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From  equations  (7),  (8),  (10)  and  (13): 


K*  £  =  u*  =  u*  -  u 
P  s 


X  f,  *  (x  ,  £  ,  t)  x  -  £-  (x  ,  t) 
v*  _  m  -  1  v  m  *  *  J  m  2  ^  m  ' 

iC”  £  =  - 


g 


1 


g- 


(14) 


X  -  f  *  (x  ,  £  ,  t)  X  *  f,  (x  ,  t  ) 

K*  =  m  1 .  m  /  .  m  .  2  m  .  (15) 


e  Si 


§2 


if 


Equation  (15)  gives  us  the  gain  which  will  insure  that, 
xp  =  xm  1  e  initially  (Figure  4) ,  the  proper  control 


will  be  applied  to  maintain  it  there  if  the  system  is 
undisturbed.  It  is  expressed  solely  in  terms  of  the  model 
state  and  state  rate,  the  specified  error,  and  the  func¬ 
tional  relationships  which  define  the  plant  and  model 
responses . 
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FIGURE  4.  PLANT  AND  MODEL 

STATE  TRAJECTORIES 

The  use  of  K* 

gain  is  illustrated  in  the  following 

example : 

Example  1:  Linear, 

first  order,  time  invariant  model  and 

plant . 

X 

=  a  x„+  u  > 

P 

P  P 

X 

1 

=  b  x  +  u  ' 

/  a  f  b 

1 

m 

m  s  ; 

we  have  x  * 

=  a  x  *  ♦  u  * 

From  equations  (5)  and  (6): 


FIGURE  5.  CONTROL  SCHEMATIC  FOR  EXAMPLE  1 


The  system  was  simulated  on  the  IBM  3033  using  the 
Continuous  system  Modeling  Program  (CSMP) .  The  following 


values  were  used: 


a  =  +  0 . 5 


b  =  '  1.0 


e  =  Q  .  0  5 


•  0 


XptO)  =  0.1 


(unstable  impulse  response) 


DQT  {PLANT)  =  0.  5*X  (PJ, ANTI  *U  (PLAN1 


'  8t 


*b.oo  o.s 


l"  '"[ME  *•“  '■* 


FIGURE  5a.  SIMULATION  RESULTS  FOR  EXAMPLE  1 


This  value  of  e  was  chosen  so  that  the  error 
following  could  be  shown  graphically.  The  system  performs 
similarly  for  any  value  of  e  >  0  . 


The  input  (us)  was  a  unit  step  at  t  =  0  .  A  time 
history  of  the  system  response  for  the  first  five  seconds  is 
shown  in  Figure  5a. 

B.  STABILITY  OF  THE  SYSTEM  USING  K*e  CONTROL 

The  question  to  be  answered  is:  Will  the  gain  K* 

satisfactorily  control  the  plant  if  e  ^  e  ? 

We  require  that,  in  response  to  an  error  such  that 

e  f  e  ,  the  system  tend  toward  e  =  e  .  That  is,  the  state 

of  the  plant,  x  ,  should  tend  toward  x  +  e  (static 

p  m  — 

stability)  . 


FIGURE  6.  PLANT  TRAJECTORIES  FOR  DECREASING  THE  ERROR 


If  e  >  e  ,  we  require  6  <  0  : 

e  >  (xB  •  ip)  sgn  (xm  -  xp)  <  0  if  e  >  e  (1?) 

(Again,  the  Signum  function  has  no  derivative  since  the 

plant  state  trajectory  does  not  cross  that  of  the  model.) 

We  have  x  *  =  x  (6)  so, 
p  m  v  J  ’ 


6  =  (xp*  -  5cp)  sgn  (xm  -  xp)  <  0  if  e  >  £ 


(18) 


The  two  cases  are  x  >  x  and  x  <  x 

m  p  m  p 


X  >  X  =>  X  *  <  X 
m  p  p  p 


for  6  <  £ 


(19) 


X  X  =>  x  *  x„ 
m  <  p  p  >  p 

These  two  cases  are  shown  in  Figure  6. 
From  equation  (1) : 


*p  '  fl  (xp  .  t)  *  S!  K«  e  | 


*  = 


fl  (Xn*  »  t)  +  Si  K< 


(20) 


In  determining  whether  x  *  <  x  or  k  *  >  k  ,we 

P  P  P  P 

know  e  >  e  ,  and  we  can  examine  f^  (x  ,  t)  to  see  if 
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it  increases  or  decreases  with  xp  .  But  the  sign  and 
magnitude  of  K*  will  be  determined  by  the  model  being 
followed  (equation  (15)). 

From  equation  (20)  we  have 

xp  *  *p*  -  [fx  Up  ,  t)  -  fx  (xp*  ,  t)  ]  +  gx  K*  (e  -  e)  (21) 

From  equation  (19)  ,  with  e  >  e  , 

First  Case:  x_  >  x_  =>  x„*  <  x„  => 

-  m  p  p  p 

[f1  (Xp  ,  t)  -  (xp*  »  t)  ]  +  gx  K*  (e  -  e)  >  0 

or  [f:  (xp  ,  t)  -  f2  (xp*  ,  t) ]  >  ■  gj  K*  (e  -  e)  (22a) 

Likewise, 

Second  Case:  X  <  x„  => 

-  m  p 

[fl  (xp  ,  t)  -  fx  (Xp*  ,  t)]  <  -  gl  K*  (e  -  e)  (22b) 

For  any  given  xp  ,  xm  ,  e  and  g-^  ,  all  the  quantities 

in  equations  (22a)  and  (22b)  except  K*  are  determined. 
Since  we  wish  our  descriptions  of  the  plant  and  the  model 
to  be  arbitrary,  we  cannot  assure  that  equations  (22a)  and 
(22b)  will  hold  for  all  cases. 
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This  is  illustrated  by  Example  1,  where 


x  =  a  x„  +  u 
P  P  P 


x  =  b  +  u 
m  ms 


Using  equations  (22a)  and  (22b) 


x  >  x  requires 

m  p 


For  x  >  x  , 
m  p 

since  (e  -  e) 


a  xp  -  a  Xp*  >  -  K*  (e  -  e) 
a  (Xp  -  Xp*)  >  -  K*  (e  -  e) 
(xp  -  Xp*)  =  -  (e  -  e) 

>  0  ,  we  require 

-  a  >  -  K*  or  a  <  K* 


Similarly,  for  x  <  x  we  require 

m  p 

a  <  -  K* 

Thus,  if  a  <  0  ,  we  require  | a j  <  |K*|  ,  or 

|a|  <  |  (b  '  a)  xa  *  a  e  *6"  (*„  ~  *p*> 

using  equation  (15) .  This  places  an  unacceptable 
tion  on  the  choice  of  model  and  e  .  If  a  >  0  , 
the  other  of  the  inequalities  is  not  satisfied. 


restric- 
one  or 
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Using  the  values  given  in  Example  1 ,  the  value  of  K* 
at  t  =  1.5  seconds  was  found  (from  the  simulation)  to  be 
-  12.2  .  Here,  a  =  0.5  ,  a  is  not  less  than  K*  and 
we  can  expect  divergence  if  xm  >  xp  • 

A  negative  step  disturbance  of  magnitude  0.3  was  super 
imposed  on  Xp  at  t  =  1.5  seconds  in  the  simulation.  T 
resulting  divergence  is  shown  in  Figure  7. 


FIGURE  7.  DIVERGENCE  RESULTING  FROM  K*e  CONTROL 


C.  STABILIZING  THE  SYSTEM 

Neutral  static  stability  (6  =  0)  can  be  established  by 
taking  e  =  e.  That  is,  if  an  error  (e)  is  present  in 
the  system,  we  take  that  error  as  our  new  e  .  If  the 
error  changes,  we  change  e  .  Equation  (IS)  becomes 


kn  ■ 


Sn  '  £1N  (xm  '  e  >  tJ  *m  ’  f2  (xm  • 


e  g- 


e  §2 


(23) 


where  the  subscript  N  denotes  a  neutrally  stable  system, 
and 


^1N  ^xm  *  e  »  >  t) 


1  v  pN 


x  ..  =  x  -  e  sgn  (x  -  x  ) 
pN  m  s  m  p 


Since  uKN  =  e 


u 


KN 


.  *m  '  f IN  txm  ■  e  ■  tJ 

8, 


c  -  f7  (x  ,  t) 
m  2  v  m  * 

g2  ' 


(24) 


At  this  point  we  are  no  longer  computing  a  gain,  but 
directly  synthesizing  a  control  modification  using 
xm  ,  xm  ,  e  and  the  descriptions  of  plant  and  model 
dynamics . 

Equation  (24)  gives  the  control  which,  when  added  to 
the  system  control,  will  cause  the  plant  to  follow  the 
model  with  constant  error.  If  the  error  present  in  the 
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system  should  change,  the  plant  will  follow  the  model 
with  the  new  error  held  constant.  This  is  illustrated  by 
way  of  Example  1.  Equation  (16)  becomes 


FIGURE  8.  CONSTANT  ERROR  CONTROL  SIMULATION 
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The  simulation  was  run  with  x  being  disturbed  bv 

p 

step  inputs  at  t  =  1.0  and  t  =  2.0  seconds  as  shown 
in  Figure  8. 

Note  that  the  error  is  constant  even  with  e  =  0  . 

This  does  not  alter  the  discussion  following  Figure  1 
regarding  steady  state  zero  error,  because  we  are  no 
longer  calculating  a  gain  as  shown  in  Figure  1.  The  present 
form  of  the  system  is  shown  in  Figure  9. 


FIGURE  9.  FORM  OF  THE  SYSTEM  WITH  CONTROL  SYNTHESIS 
At  this  point,  if 

A.  The  physical  pant  is  perfectly  described  by  the 
plant  equation,  and 

B.  The  plant  and  the  model  have  the  same  initial 
conditions  (e(0)  *  0)  ,  and 


C.  There  are  no  external  disturbances  (or  if  the  mean 
disturbance  is  zero) ,  then  the  control  defined  by  equation 
(24)  will  yield  zero  error  (or  zero  mean  error) . 

D.  ERROR  AND  ERROR  RATE  CONTROL 

If  an  error  exists  in  the  system,  it  is  desired  that 
it  be  reduced  to  zero  in  a  controlled  manner.  Figure  10 
shows  a  system  with  initial  error.  The  trajectory  of  Xp 
with  constant  error  (equation  (24))  is  shown,  as  well  as 
a  desired  response  to  this  condition  (dashed  line) . 


FIGURE  10.  CONSTANT  ERROR  AND  DESIRED  RESPONSE 
TRAJECTORIES  OF  THE  PLANT 

The  desired  response  may  be  obtained  by  observing  that 
the  plant  may  be  made  to  follow  any  model  with  constant 
error.  That  is,  if  we  can  find  another  model  (to  be  called 
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the  control  model)  which,  if  followed  with  constant  error, 
will  cause  the  plant  state  to  return  to  the  model 
trajectory,  the  desired  response  is  obtained. 

The  trajectory  of  the  control  model  is  shown  in 
Figure  11. 


FIGURE  11.  CONTROL  MODEL  TRAJECTORY 


The  notations  e  and  e  represent  the  errors 

m  cm  r 

between  the  plant  and  the  model,  and  between  the  plant  and 

the  control  model,  respectively.  Note  that  ecm  is  held 

constant,  and  that  the  trajectory  of  x  causes  e„  to 

J  cm  m 

diminish . 

From  Figure  10  the  nature  of  x_  is  seen:  If  e  is 
&  cm  m 

0  ,  then  x  is  the  same  as  x  (except,  perhaps,  for 
cm  m 

initial  conditions).  If  e  is  non-zero,  then  x _ 

m  cm 

differs  from  xm  by  an  amount  proportional  to  em  . 


A  control  model  which  achieves  the  desired  response  may 
be  written: 

*cm  ‘  *m  4  (V>  s«n  (*m  '  ltp)  l26> 

where  A  (e„)  >  0  and 

m  — 

A  (em)  *  0  if 


e  =  0 
m 


The  Signum  function  assures  that 


x„„  <  x  if  x  >  x  and 
cm  m  pm 


x_„  >  *  if  x  <  xm,  as  required, 

cm  m  pm  ^ 


From  equation  (2) : 


:cm  '  f2  (xm  •  tJ  *  «2  us  *  4  temJ  sgn  (  xm  ‘  V  (27) 


so 


u_ 


cm 


-  f7  C\  ,  t)  -  A  (ej  sgn  (x, 


m 


m 


xp) 


§7 


(28) 
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The  additive  control  input  necessary  to  follow  the 


control  model  with  constant  error  is  then  (from 
equation  (24)): 


xcm  '  flN  ^xcm  *  ecm  ’  t) 

UK  ^ 


(29) 


f,  (xm  ,  t)  -  A  (em)  sgn  (xm  -  xp) 


Lcm  2  1  m 


In  equation  (29)  f,.r  (x  .  e„  ,  t)  is  defined  as 
n  IN  cm  cm 

the  function  evaluated  at 


xp  ■  xcm  '  ecm  SS"  lxcm  '  xp> 


This  expression  may  be  shown  to  be  independent  of  xcm 
by  observing  the  six  possible  relationships  between 
x  ,  x  and  x  and  writing  the  expression  for  e 


p  m 
for  each 


cm 


cm 


x„  >  xm  >  x„„ 
p  m  cm 


e  =  x  x _ +  e 

cm  m  cm  m 


X„  >  x, _ >  xm 

p  cm  m 


e  *  x  x _ +  e 

cm  m  cm  m 


X  >  X  >  X 

m  p  cm 


e_  =  x  ~  e  -  e 
cm  m  cm  m 


x„  >  x^  >  X, 
m  cm  i 


e  *  -  x  +  x  +  e 
cm  m  cm  m 


X  >  x„  >  xm 
cm  p  m 


e  =  -  x  +  x  -  e 
cm  m  cm  m 


>  x„  >  xm 
cm  m  p 


cm 


‘  x„  +  x„  +  e„ 
m  cm  m 
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•'1 


From  which,  for  all  cases 


e; _ =  x„m  sgn  (x.„  -  x  )  -  x  sgn  (x  -  x  ) 

cm  cm  3  cm  p  m  6  cm  p 

+  e  sgn  (xm  -  x  1  sgn  (x, _ -  x  ) 

m  ®  m  p  &  cm  p' 

then 

e, _ sgn  (x, _ -  xl  =  x, _ -  x  +  e  sgn  (x  -  x  ) 

cm  6  cm  p'  cm  m  m  3  v  m  p 

so  x, _ -  e _ sgn  (x, _ -  x  )  =  x  -  e  sgn  (x  -  x  ) 

cm  cm  &  v  cm  p'  m  m  3  v  m  p' 

The  expression  for  is  then 


§2 


where  f,v  (x  ,  e  ,  t)  is  the  function  f,  evaluated 
i  m  m  1 

at  x  =  x  -  e  sgn  (x  -  x  )  . 

p  m  m  6  v  m  p 

In  equation  (30)  ,  if  g^  =  g?  then  u^  does  not 

depend  on  x 
r  cm 

In  any  case,  it  is  simpler  at  this  point  to  replace 
the  last  term  in  equation  (30)  by  ug  and  calculate 
directly  u  *  u„  +  u  ,  or 

p  K  S 
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(31) 


u 


P 


*cm  f IN  ^xm  » 

Si 


e 

m 


t) 


The  selection  of  the  function  A  (e  )  appears  to  be 

arbitrary  with  unconstrained  control.  The  plant  following 

the  control  model  in  exactly  the  desired  manner  in  each 

simulation  run.  The  form  of  A  (e  )  chosen  for  most 

v  nr 

simulations  was  A  (e  )  =  [exp  (w  .  e  )  -  1]  ,  where  w 

m  r  m 

was  selected  to  vary  the  speed  of  the  response. 

The  application  of  this  form  of  error  control  is 
illustrated  by  continuing  Example  1 : 

a  f  b 

From  equation  (26),  *cm  =  b  cm  +  ug  +  A  (ej  sgn  (xm  -  x  ) 


From  equation  (30) ,  uR 


cm 


a  (xm  -  e  sgn  (x, 


m 


m 


m 


cm 


-  b  x. 


m 


-  4  (en)  sgn  (xm  -  x  ) 

- r - - 
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simplifies  to 


The  expression  for 


u 


K 


b  x  + 


m 


(e  ) 
v  nr 


(x 


m 


•  a  lxm  '  em  s*n  (xm  ‘  xp>> 

selecting  A  (e  )  =  [exp  (3  e )  -  1]  . 
o  ^  m  r  m 

The  problem  was  simulated  with  the  values  a  *  +  0.5  , 
and  b  *  -  1.0  as  before.  The  plant,  model,  and  control 
model  had  the  same  initial  conditions  (zero).  A  step 
disturbance  was  imposed  on  the  plant  at  t  =  0.5  second. 
The  response  is  shown  in  Figure  12: 


FIGURE  12.  SIMULATION  USING  CONSTANT  ERROR 
MODEL  FOLLOWING  CONTROL 
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The  system  at  this  point  is  most  simply  described 


from  equation  (31)  as  a  feedforward  control  synthesizer  with 
error  feedback,  as  shown  in  Figure  15: 


FIGURE  13.  CONTROL  SYSTEM  CONFIGURATION  WITH 

CONSTANT  ERROR  MODEL  FOLLOWING  CONTROL 


■i 

1 


V.  SUMMARY 


The  method  of  constant  error  model  following  control 
is  summarized  as  follows: 

Given  the  plant  described  by 

*p  =  fl  (xp  »  +  §1  up  (1) 

and  the  model 


*m  *  h  txm  •  tJ  *  *»  u 


2  s 


(2) 


Define  the  control  model 


>c  «  x  +  A  (em)  sgn  U  -  x  ) 
cm  m  m  &  m  p 


(26) 


Select  A  (e  )  such  that 
m 


A  (e  )  >  0 
m  — 


A  (0)  *  0 


A  (eml)  >  A  (em7)  if  em1  >  e 


m2J 


'ml  m2 


Solve  equation  (1)  for  Up  : 


u 


*P  -  h  (XP  ■  1  } 

Si 


(32) 
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In  equation  (32)  substitute 


x 


cm 


fl  (xp  »  t;)  =  £1N  (xm  ’  em  »  ^ 


where 


■IN 


Cxm  • 


e  ,  t) 
m  *  J 


(x  ,  t)  with 


x 


P 


e  sgn  fx 
m  °  ^  m 


The  result 


UP  = 


x_m  -  f1v  (xm  ,  em  ,  t) 
cm  IN  m  m 


(31) 


Is  the  plant  input  which  makes  the  plant  states  equal  to 
those  of  the  model,  and  restores  this  condition  in  the 
presence  of  error. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  It  is  concluded  that  constant  error  model  following 
control  achieves  the  desired  goal  of  error  control  for  the 
systems  described  in  this  paper.  It  offers  the  advantage 
of  perceptual  simplicity  in  that  the  designer  can  visualize 
exactly  the  effect  of  his  control  method  on  the  system 
response.  It  affords  flexibility  in  that  the  choice  of  the 
restoring  function  A  (e  )  is  arbitrary  within  the  few 
constraints  mentioned. 

B.  The  following  areas  of  future  research  are  suggested: 

1.  Development  of  the  theory  with  application  to  highe 
order  systems  with  non-additive  controls. 

2.  Investigation  of  the  response  of  systems  in  which 
the  physical  plant  is  not  accurately  described  by  the  plant 
equations . 

3.  Application  of  the  theory  to  practical  problems. 


'APPENDIX  A:  EXAMPLE  OF  APPLICATION  OF  THE  METHOD  TO 
A  FIRST  ORDER  NONLINEAR  SYSTEM  WITH 
TIME  VARYING  COEFFICIENTS 


1.  Assume  the  Plant  is  given  by 


x 


P 


+ 


2up 


(Al) 


and  the  model  by 


x  = 
m 


(t  x  )1/2  ♦  u 
m  s 


(A2) 


define 


cm 


=  *  +  A 

m 


s«n  (xm  '  V 


(AJ) 


then 


u 


P 


into  (AS)  substitute 


(A4) 


(A5) 


xp  ‘ 


m 


-  e_  sgn  (xm  -  xp) 


m 


(A6) 
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then 


(A7) 

2.  The  system  was  simulated  using  the  Continuous  System 

Modeling  Program  (CSMP)  on  the  IBM- 3033.  Time  histories 
of  the  response  following  a  disturbance  at  t  =  0.5 
second  are  shown  in  Figures  A1  through  A8 .  The  system 
control  input  is  a  unit  step  at  t  =  0  in  all  cases. 
Shown  also  are  the  plant  controls  generated  in  each  case. 
The  four  cases  were  the  same  except  for  the  choice  of 

ACem]  : 

A.  A (em)  *  5em  (Figures  A1  and  A2) 

A(em)  ■  exp  (em)  -  1  (Figures  A3  and  A4) 

A (em)  *  exp  (2em)  "  1  (Figures  A5  and  A6) 

A(ejn)  =  exp  (3era)  -  1  (Figures  A7  and  A8) 
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FIGURE  Al.  MODEL  AND  PLANT  RESPONSES, 
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FIGURE  A7 :  MODEL  AND  PLANT  RESPONSES,  A ( e  )  *  exp  (3  e  ) - 1  . 
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